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ABSTRACT 

We study the amount and distribution of dark matter substructures within dark matter 
haloes, using a large set of high-resolution simulations ranging from group size to 
cluster size haloes, and carried our within a cosmological model consistent with WMAP 
7-year data. In particular, we study how the measured properties of subhaloes vary 
as a function of the parent halo mass, the physical properties of the parent halo, and 
redshift. The fraction of halo mass in substructures increases with increasing mass: 
it is of the order of 5 per cent for haloes with M200 ~ lO^'^ M© and of the order 
of 10 per cent for the most massive haloes in our sample, with M200 ^ 10^^ M©. 
There is, however, a very large halo-to-halo scatter that can be explained only in 
part by a range of halo physical properties, e.g. concentration. At given halo mass, 
less concentrated haloes contain significantly larger fractions of mass in substructures 
because of the reduced strength of tidal disruption. Most of the substructure mass is 
located at the outskirts of the parent haloes, in relatively few massive subhaloes. This 
mass segregation appears to become stronger at increasing redshift, and should reflect 
into a more significant mass segregation of the galaxy population at different cosmic 
epochs. When haloes are accreted onto larger structures, their mass is significantly 
reduced by tidal stripping. Haloes that are more massive at the time of accretion 
(these should host more luminous galaxies) are brought closer to the centre on shorter 
time-scales by dynamical friction, and therefore suffer of a more significant stripping. 
The halo merger rate depends strongly on the environment with substructure in more 
massive haloes suffering more important mergers than their counterparts residing in 
less massive systems. This should translate into a different morphological mix for 
haloes of different mass. 

Key words: cosmology: dark matter - clusters: general - galaxies: evolution - galaxy: 
formation. 



1 INTRODUCTION 

In the currently accepted ACDM paradigm for cosmic struc- 
ture formation, small dark matter haloes form first while 
more massive haloes form later through accretion of dif- 
fuse matter and mergers between smaller systems. During 
the last decades, we have witnessed a rapid development of 
numerical algorithms and a significant increase in numeri- 
cal resolution, that have allowed us to improve our knowl- 
edge of the formation and evolution of dark matter struc- 
tures. In particular, the increase in numerical resolution 
has allowed us to overcome the so-called overmerging prob- 
lem, i.e. the rapid disr uption of galaxy-siz e substructures 
in groups and clusters (|Klvpin et al.l 1 19991 . and references 
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therein). If any, we are now facing the opposite problem, 
at least on galaxy scales, where many more substructures 
than visible dwarf g alaxie s are found (Ishivama ct al.. . 20091 : 
Prikhonov fc Klvpinll2009l . and references therein). 

According to the two stage theory proposed by 
IWhite fc Reed l|l978h , the physical properties of galaxies are 
determined by cooling and condensation of gas within the 
potential wells of dark matter haloes. Therefore, substruc- 
tures represent the birth-sites of luminous galaxies, and the 
analysis of their mass and spatial distribution, as well as of 
their merger and mass accretion histories provide important 
information about the expected properties of galaxies in the 
framework of hierarchical galaxy formation models. 

Nowadays, a wealth of substructures are routinely iden- 
tified in dissipationless simulations, and their statistical 
properties and evolution have been studied in detail in the 
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past years. The identification of dark matter substructures, 
or subhaloes, remains a difficult technical task that can be 
achie ved using different algorithms (see e.g. iKnebe et al.l 
Each of these has its own advantages and weaknesses, 
and different criteria for defining the boundaries and mem- 
bership of substructures are likely leading to systematic dif- 
ferences between the physical properties of subhaloes iden- 
tified through different algorithms. However, these might be 
probably corrected using simple scaling factors, as suggested 
by the fact that different studies find very similar slopes 
for the subhalo mass function, i.e. the distribution of sub- 
structures as a function of their mass. This is one of the 
most accurately studied properties of dark matter substruc- 
tures, although it rema ins unclear if and h ow it depends on 
the parent halo mass. iMoore et al.l l| 19991 ) used one high- 
resolution simulation of a cluster-size halo and one high- 
resolution simulation of a galaxy-size halo, and found that 
the latter can be viewed as a scaled version of the former. 
Later work bv lDe Lucia et akl l|2004l 'l used larger samples of 
simulated haloes, but found no clear variation of the sub- 
halo mass function as a function of the parent halo m ass. 
Such a dependen cy was later found by Gao ct al] (120041 ) and 
iGao et al.l HoTl), who showed that the subhalo mass func- 
tion varies systematically as a function of halo mass and halo 
physical properties like concentration and formation time. 

Typically, only about 10 per cent of the total mass of a 
dark matter halo is found in substructures. In addition, their 
spatial distribution is f ound to be anti-bia s ed with respect 
to that of dark matter (iGhigna et al.ll2000l: iDe Lucia et al.l 
l2004l : lNagai fc Kravtsov 2005; S aro et al.ll2010l ). It is unclear 
if the radial distribution of subst r uctur es depends on the 
parent halo mass. iDe Lucia et al.l (|2004l ) found hints for a 
steeper radial number density profiles of substructures in low 
mass haloes than in high mass haloes. They used, however, 
a relatively small sample of simulated haloes, that were run 
with different codes and numerical parameters. In this study, 
we will re-address this issue by using a much larger sample of 
simulated haloes, all run with the same code and numerical 
parameters. 

Most previous work focusing on dark matter substruc- 
tures has studied their properties as a function of their 
present day mass. This quantity cannot be, however, sim- 
ply related to the luminosity of the galaxies residing in the 
substructures under consideration. Indeed, dark matter sub- 
structures are very fragile systems that are strongly affecte d 
by tidal stripping l|De Lucia et al.l |2004| : iGao et all [2004 ) . 
Since this process affects primarily the outer regions of sub- 
haloes, and galaxies reside in their inner regions, it is to 
be expected that the galaxy luminosity/stellar mass is more 
strongly related to the mass of the substructure at the time 
of infall (i.e. befo r e becoming a substructure) than at presen t 
l|Gao et al.ll2003 : IWang et al.ll2006l : I Vale fc Ostrikeifl2006l ). 
In this paper, we will study the evolution of dark matter sub- 
structures splitting our samples according to different values 
of the mass at infall. 

In this paper, we take advantage of a large set of N- 
body simulations covering a wide dynamical range in halo 
mass, and with relatively high-resolution. This will allow 
us to study how the statistical properties of substructures 
vary as a function of halo mass, cosmic epoch, and physical 
properties of the parent halo. The layout of the paper is as 
follows: in section [2] we introduce the simulation set and 



samples used in our study. In section [S] we study how the 
subhalo mass function and subhalo spatial distribution vary 
as a function of halo mass, redshift and concentration. In 
the second part of our paper (section|4ll, we discuss the mass 
accretion and merging histories of subhaloes as a function 
of their mass, accretion time, and environment. Finally, in 
Section O we discuss our findings and give our conclusions. 



2 CLUSTER SIMULATIONS 

Our set of DM haloes is based on 'zoom-in' simulations of 
27 Lagrangian regions extracted around massive dark matter 
haloes, originally identified within a low-resolution N-body 
cosmological simulation. F or a detailed di s cussio n of this 
simulation set, we refer to lBonafed"e^ et all (|201ll . see also 
iFabian et aLll201ll ). The parent simulation followed 1024^ 
DM particles within of a box of l/i~^Gpc comoving on a 
side. The adopted cosmological model assumed Qm = 0.24 
for the matter density parameter, fibar = 0.04 for the con- 
tribution of baryons. Ho = 72kms~^Mpc~^ for the present- 
day Hubble constant, Us = 0.96 for the primordial spec- 
tral index, and as = 0.8 for the normalization of the power 
spectrum. The latter is expressed as the r.m.s. fluctuation 
level at z = within a top-hat sphere of 8/i~'^Mpc ra- 
dius. With this parameters choice, the assumed cosmogony 
is consistent with constraints derived from seven-year data 
from the Wilkinson M icrowave Anisotropy Probe (WMAP7, 
iKomatsu et al.ll201ll ). 

The selected Lagrangian regions were chosen so that 
13 of them are centred around the 13 most massive 
clusters found in the cosmological volume, all having 
viriaQ mass M200 — 10^^ h~^MQ. Additional regions 
were chosen around clusters in the mass range M200 — 
(5 - 10) X 10"/i"^A4"q. Within each Lagran gian region, 
we increased mass resolution and added the relevant 
high-frequency modes of the power spectrum, using the 
Zoomed Initial Condition (ZIC) technique presented by 
iTormen. Bouchet fc White! l| 19971 ). Outside the regions of 
high-resolution, particles of mass increasing with distance 
are used, so that the computational effort is concentrated 
on the cluster of interest, while a correct description of 
the large-scale tidal field is preserved. For the simulations 
used in this study, the initial conditions have been gener- 
ated using TTiDM ~ 10* h~^MQ for DM particle mass in the 
high-resolution regions. This m ass resolution is a fact or 10 
better than the va lue used by iBonafede et al.l l|201ll ) and 
iFabi an et all (|201ll ) to carry out hydrodynamic simulations 
for the same set of haloes. 

Using an iterative procedure, we have shaped each high- 
resolution Lagrangian region so that no low-resolution par- 
ticle 'contaminates' the central 'zoomed in' halo, out to 5 
virial radii of the main cluster at z = 0. In our simulations, 
each high resolution region is sufficiently large to contain 
more than one interesting massive halo, with no 'contami- 
nants', out to at least one virial radius. Our final sample con- 
tains 341 haloes with mass larger than lQ^''^h~^ Mq. We have 
split this sample into 5 different subsamples, as indicated 

^ Here we define the virial mass (M200) as the mass contained 
within the radius -R200 1 that encloses a mean density of 200 times 
the critical density of the Universe at the redshift of interest. 
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Table 1. Our simulation set has been split in five subsamples, 
according to the halo mass. In the first column, we give the name 
of the subsample, while the second column indicates the range 
of M200 values corresponding to each subsample. The third and 
fourth columns give the number of haloes and mean number of 
subhaloes (with mass above 2 ■ 10^ /i~^Mq) within the virial ra- 
dius {-R20o)i respectively. 



Name 


Mass range 


^haloes 


* subs 


SI 




13 


2943 


S2 


[5-10] X lO^^h-'^MQ 


15 


1693 


S3 


[1-5] X IQi"* /i-1Mq 


25 


358 


S4 


[5-10] X IQis h-iMg 


29 


146 


S5 


[1-5] X 10^3 /i-1Mq 


259 


40 



in Table [T] where we list the number of non-contaminated 
haloes for each sample and the mean number of substruc- 
tures per halo in each subsample. 

Simulations have been carried out using the Tree-PM 
GADGET-3 code. We adopted a Plummer-equivalent soft- 
ening length for the computation of the gravitational force in 
the high-resolution region. This is fixed to e = 2.3 kpc 
in physical units at redshift z < 2, and in comoving units at 
higher redshift. For each simulation, data have been stored 
at 93 output times between z ~ 60 and z — 0. Dark mat- 
ter haloes have been identified using a standard friends- 
of-friends (FOF) algorithm, with a linking length of 0.16 
in units of the mean inter-particle separatio n in the high- 
resolu tion region. The algorithm SUBFIND (jSpringel et al.l 
I2OOII ) has then been used to decompose each FOF group 
into a set of disjoint substructures, identified as locally over- 
dense regions in the density field of the background halo. As 
in previous work, only substructures which retain at least 
20 bound particles after a gravitational unbinding proce- 
dure are considered to be genuine substructures. Given our 
numerical resolution, the smallest structure we can resolve 
has a mass oi M — 2x 10® h~^MQ. To avoid being too close 
to the resolution limit of the simulations, we will sometimes 
consider only substructures that contain at least 100 parti- 
cles, i.e. we will adopt a mass limit of 1 x 10^'' h~^MQ. 



3 AMOUNT AND DISTRIBUTIONS OF DARK 
MATTER SUBSTRUCTURES 

In this section we will consider some basic statistics of the 
dark matter substructures in our sample. In particular, we 
will address the following questions: what is the mass frac- 
tion in substructures? What is their mass and spatial dis- 
tribution? And how do these properties vary as a function 
of the halo mass, or as a function of other physical proper- 
ties of the parent haloes? As discussed above, if subhaloes 
are to be considered the places where galaxies are located, 
these statistics provide us important information about the 
statistical properties of cluster galaxy populations expected 
in hierarchical cosmologies. 

3.1 Mass Fraction in Subhaloes 

Previous work has found that only 5 to 10 per cent of 
the halo mass is contained in substructures, with most of 



it actually co ntained in relatively few mass i ve subhaloes 
jGhigna et al.|[l 998 200 ol: ISpringel et al.ll200ll : IStoehr etHI 
l2003l : lGao eral.ll2004l : IDe Lucia et al.ll2004D . 



Results for our simulation set are shown in Figure [T] 
The top left panel shows the cumulative mass fraction in 
subhaloes above the mass indicated on the x-axis, for the 
five samples considered in our study. There is a clear trend 
for an increasing mass in substructures for more massive 
haloes. For our most massive sample (SI), about ten per 
cent of the halo mass is contained in substructures more 
massive than 2 x 10^ Mq, and approximately ten per 
cent of the mass in substructures is contained in the most 
massive ones. For less massive haloes, the mass fraction in 
substructures decreases. 

Most of the substructures are located outside the cen- 
tral core of dark matter haloes. In particular, the top right 
panel of Figure [T] shows that the substructure mass fraction 
is smaller than ~ 1 per cent out to ~ 0.3 x r2oo , and increases 
to half its total (within r^oo) value at ~ 0.8 x r2oo- The re- 
sults shown can be explained by considering that haloes of 
larger mass are less concentrated and dynamically younger 
than their less massive counterparts. As we will show below, 
and as discussed in previous studies, subhaloes are strongly 
affected by dynamical friction and tidal stripping. Less mas- 
sive haloes assemble earlier than their more massive coun- 
terparts, i.e. accrete most of the haloes that contribute to 
their final mass at early times, so that there was enough 
time to 'erase' the structures below the resolution of the 
simulation in these systems. In addition, haloes that were 
accreted earlier, and therefore suffered of tidal stripping for 
longer times, are preferentially lo cated closer to the centre 
(see Figure 15 in lGao et alll2004l '). 

For haloes of the same mass, a relatively large range of 
concentrations is possible so that a range of mass fractions is 
expected. This is confirmed in the bottom panels of Figure[T] 
where we have considered only haloes in our least massive 
sample (S5), and split it into three different bins according 
to the halo concentration so as to have the same number of 
haloes for each bin. We approximate the concentration by 
Vmax/V2oo, where Vmax is the maximum circular velocity, 
which is computed by cons idering all pa rticles bound to a 
given halo, while V200 ~ a/ GM200 / -R200 • Interestingly, the 
lowest concentration bin contains substructure mass frac- 
tions that are, on average, very close to those of our most 
massive samples (SI in the top panels). This confirms that 
the halo to halo scatter is very large, and that it can be 
explained only in part by haloes in the same mass bin cov- 
ering a range of physical properties. In order to give an idea 
of the intrinsic scatter of haloes in the same mass bin, we 
have repeated the last point in the top right panel of Fig- 
ure [TJ showing this time the median and the 25th and 75th 
percentile of the distributions obtained at R/R200 = 1- 

3.2 Subhalo Mass Function 

One of the most basic statistics of the subhalo population is 
provided by the subhalo mass function, i.e. the distribution 
of dark matter substructures as a function of their mass. 
This has been analysed in many previous studies with the 
aim to answer the following questions: does the subhalo mass 
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Figure 1. Top panels: cumulative mass fraction in substructures as a function of subhalo mass (left) and normalized distance from the 
halo centre (right), for the five samples used in this study (different symbols, as indicated in the legend). In the right panel the rightmost 
symbols with error bars show the median, 25"* and percentile of the distributions at R/R200 = 1- Bottom panels: same as in the 

top panels but using only haloes from our sample S5 (the least massive one), and splitting the sample in three different bins according 
to the concentration of the parent haloes. In all panels, symbols connected by lines show the mean values, while error bars show the rms 
scatter around the mean. 



function vary as a function of the parent halo mass? How 
does it vary as a function of cosmic time? And as a function 
of halo properties (e.g. concentration, formation time, etc.)? 

First studies were based on very small samples of sim- 
ulated haloes, and claimed the 'univers ality' of the subhalo 
mass function. E.g. iMoore et al.l ()l999l ) compared the sub- 
structure mass distribution obtained for one simulated clus- 
ter of mass similar to that of the Virgo Cluster, and one 
simulated galaxy-size halo, and argued that galactic haloes 
can be considered as 'scaled versions' of cluster-size haloes. 
iDe Lucia et~al] l|2004) used a sample of ~ 11 high resolution 
resimulations of galaxy clusters together with a simulation 
of a region with average density. They argued that the sub- 



halo mass function depends at most weakly on the parent 
halo mass and that the (nearly) invariance of the subhalo 
mass function could lie in the physical nature of the dynam- 
ical balance between two opposite effects: the destruction of 
substructures due to dynamical friction and tidal stripping 
on the one hand, and the accretion of ne w substruct u res on 
the other hand. Contemporary work by iGao et al.l (|2004| ) 
and later work (e.g. Gao et al. 20 111) has demonstrated that 
the subhalo mass function does depend on the parent halo 
mass, as well as on the physical properties of the parent halo, 
in p articular its concen tration and formation time. We note 
that ICao et al.l (|2004l ) used a sample of simulated haloes 
that was not homogeneous in terms of resolution (typically 
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Figure 2. Slope of the differential mass function measured for 
the different samples considered in this study, at different cosmic 
epochs (solid Une for 2 = 0, dotted for z = 0.5, dashed line for 
2 = 1, and dash-dotted line for 2 = 2). Error bars are computed 
as the standard deviation of the slopes measured for each halo 
within the sample. For reasons of clarity, a small shift has been 
added to the abscissa. 



lower than ours), cosmoloRi cal parame t ers, a nd simulation 
codes. The sample used in iGao et al] l|201lh was instead 
based on a ho homogeneous set of cosmological parame- 
ters (consistent with WMAP first-year results) and included 
simulations with resolution higher than that of our sample. 
Their sample, however, did not include very massive haloes 
(~ lO'^^ Mq). It is therefore interesting to re-address the 
questions listed above using our simulation sample. 



In Figure [2j we plot the slope of the differential mass 
function obtained by fitting a power law to the mass func- 
tio ns of each sub-sample considered in our study. Follow- 
ing iDe^Tuda^^eT^l] (|2004l ). we have restricted the fit by dis- 
carding the most massive (and rarest) substructures (those 
with mass above 10^^ h~^MQ for the samples SI and S2, 
and with mass above 10^^'^ Mq for the samples S3, S4 
and S5). We find that, albeit weakly, the slope of the sub- 
halo mass function depends on the parent halo mass, and 
that there is a weak trend for shallower slopes with increas- 
ing lookback times. The best fit values we measure vary 
in the range between ~ —0.65 and ~ —0.8, in agreement 
with results fr om previous studies (e.g. iGhigna et al.|[200ol : 
iDe Lucia eTal . 2004; Gao et al. 2004). When including the 
most massive substructures in the fit, we obtain steeper 
slopes, ranging from ~ —0.91 and ~ —0.86 at redshift 2 = 0, 
but the trends shown in Figure [H a re not altered significantly. 

As explained bv iGao et al.l (1201 ll ). the dependence of 
the subhalo mass function on halo mass is a consequence 
of the fact that more massive haloes are on average less 
concentrated and dynamically younger than their less mas- 
sive counterparts. Since the strength of tidal disruption de- 
pends on halo concentration, and since haloes of a given 



mass are on average less concentrated at higher redshift, we 
also expect that the subhalo mass function depends on time. 
Figure 13 sho ws the cumulativ e subhalo mass function (nor- 
malized as in iGao et al.|[200i ) at four different redshifts in 
the left panels and for different concentrations in the right 
panels (in these panels, only haloes identified at redshift 
zero have been considered). Top and bottom panels refer 
to the haloes in the mass range [1 — 3] x lO'^^ Mq and 
[1 — 5] X 10^** h~^MQ respectively. We derive the three sub- 
samples by splitting the range of concentration in order to 
have the same number of haloes in each subsamp le. Results 
shown in Figure [3] confirms previous findings by iGao et al] 
(|201lh , and extend them to larger parent haloes masses: 
haloes at higher redshift have significantly more substruc- 
tures than those of the same mass at later times. The figure 
suggests that there is a significant evolution between 2 = 
and 2 ~ 0.5, but it becomes weaker at higher redshifts. 
We note that for the highest redshift considered, the sub- 
halo mass function does not significantly differ from that 
found at 2 ~ 1, but we not e that this could be due to poor 
statistics. iGao et all HoTl) find a similar trend for haloes 
of similar mass. At any given cosmic epoch, there is a large 
halo-to-halo scatter which is due, at least in part, to internal 
properties of the parent halo like concentration, as shown in 
the right panels of Figure |3l For the ranges of mass shown 
in Fig. [S] low concentration haloes host up to an order of 
magnitude more substructures than haloes of the same mass 
but with higher concentration. The difference between the 
different concentration bins are larger (and significant) for 
the most massive substructures. 

In order to verify that the results of our analysis are 
robust against numerical resolution, we have compared the 
cumulative sub-halo mass function obtained for the set of 
simulated halos presented here to that obtained for the same 
halos simulated at a 10 times lower mass resolution. We 
find that the two distributions agree very well to each other, 
within the mass range accessible to both resolutions. This 
confirms that both our simulations and the procedure of halo 
identification are numerically converged. 



3.3 Radial Distribution of Subhaloes 



Previous studies i Ghigna et al."2000'; ' Pe Lucia et al.ll2004l : 
iNagai fc Kravtsov| [2005 : Saro ct al. 20i3) have shown that 
subhaloes are 'anti-biased' relative to the dark matter in 
the inner regions of haloes. No significant trend has been 
found as a function of the parent halo mass, with only hints 
for a steeper profiles of subhaloes in low massive haloes 
(|De Lucia et al.ll2004l ). 



The analysis of our sample of simulated halos confirms 
previous findings that dark matter subhaloes are anti-biased 
with respect to dark matter, with no dependence on parent 
halo mass. In fact, there is no physical reason t o expect such 
a trend. We note that iDe Lucia et al.l ((2003), who found 
hints for such a correlation, used a smaller sample of simu- 
lated haloes, that were carried out using different simulation 
codes and parameters. In contrast, our simulated haloes are 
all carried out using the same parameters and simulation 
code. 
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Figure 3. Cumulative mass functions (CMF) in units of rescaled subhalo mass, and multiplied by Msui,/M200 to take out the dominant 
mass dependence. Top and bottom panels are for haloes in the mass range [1 — 3] X 10^^ Mq and [1 — 5] X 10^'* h~^MQ, respectively. 
In the left panel, results are shown for different redshifts (solid line for z=0, dotted line for z=0.5, dashed line for z=l and dash-dotted 
line for z=2). In the right panel, only haloes identified at redshift zero have been considered, and they have been split in three bins, 
according to their concentration. Only subhaloes with more than 100 bound particles have been used to build these functions. 



iNaeai fc Kravtsovl (|2005l ) find that the anti-bias is much 
weaker if subhaloes are selected on the basis of the mass 
they had at the time of accretion onto their parent halo. 
We confirm their results in Figure 31 where we show the 
radial distribution of substructures in our sample SI (the 
most massive haloes in our simulation set). The top panel 
of Figure |3] shows the radial distribution of substructures 
selected on the basis of their present day mass, while in 
the middle panel the mass of the substructure at the time 
of accretion (defined eis the last time the halo was identi- 
fied as a central halo, see below) has been used. The figure 
shows that, in this case, selecting progressively more mas- 
sive substructures reduces the anti-bias between subhaloes 
and dark matter. The bottom panel of Figure [4] shows that 
the same is obtained by discarding substructures that are 



accreted recently. The two selections tend to pick up haloes 
that suffered a stronger dynamical friction (i.e. haloes that 
were more massive at the time of accretion) or that suffered 
of dynamical friction for a longer time (haloes that were ac- 
creted earlier). As a consequence, both selections tend to 
preferentially discard subhalos at larger radii, thus bringing 
the radial distribution of subhaloes closer to that measured 
for dark matter. 

As shown above (see right panels of Figure [1} , most 
of the substruc t ure m ass is located at the cluster outskirts. 
iDe Lucia et~all l|2004D showed that this distribution is de- 
pendent on the subhalo mass, with the most massive sub- 
structures being located at larger distances from the cluster 
centr e with respect t o less m assive substructures. In partic- 
ular, iDe Lucia et al.l (|2004l ) split their subhalo population 
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Figure 4. Radial distribution of dark matter substructures be- 
longing to haloes of the sample SI. In the top panel, different lines 
correspond to different thresholds in the Mg^i,/M2oo ratio, based 
on the present-day subhalo mass. In the middle panel, the sub- 
halo mass at the time of accretion has been considered, while in 
the bottom panel different lines correspond to subhalocs accreted 
at different times. 
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Figure 5. Cumulative radial distributions for subhaloes with 
Msub/M200 > 0.01 (soUd Une) and Msub/M200 < 0.001 (dot- 
ted line) from all samples, at different redshifts. On the y-axis, 
we plot the total mass in subhaloes within a given distance from 
the centre, normalized to the total mass in subhaloes within i?200i 
for each subhaloes population. 



in two subsamples by choosing a rather arbitrary mass ratio 
between the subhalo mass and the parent halo mass (they 
chose the value 0.01 for this ratio). Our simulations exhibit 
the same trends, but we find that this can be more or less 
'significant' depending on the particular threshold adopted 
to split the sample. In Figure [S] we show the radial dis- 
tribution of substructures with Msub/M2oo > 0.01 (solid 
lines) and Ms^i/K'hoo < 0.001 ( dashed lines). Our tre nds 
are weaker than those found by iDe Lucia et al. I l|2004 ) at 
redshift zero, when the same division is adopted. We note, 
however, that these trends are dominated by the most mas- 
sive substructure and are, therefore, significantly affected 
by low number statistics. Figure [5] also shows that the mass 
segregation becomes more important at increasing redshift. 

Considering that haloes of a given mass are less cen- 
trally concentrated and dynamically younger than their 
counterparts at later redshift, the trend found can be ex- 
plained as follows: the 'younger' haloes have massive sub- 
haloes preferentially in their outer regions because stripping 
has not had enough time to strip their outer material and 
eventually disrupt them. In more dynamically evolved clus- 
ters (those at present time), stripping has had more time 
to operate and to wash out any difference between the two 
distributions. In this picture, the balance between dynam- 
ical friction and stripping on one hand, and the accretion 
of new subhaloes on the other hand is such that the latter 
effect is dominating over the former. This is in agreement 
with the results shown above for the evolution of the cumu- 
lative mass function of substructures, whose normalization 
increases with increasing redshift. 

We stress that in Figure [5] we are considering subhaloes 
of different mass at the time they are identified. As dis- 



8 E. Contini et al. 



cussed in Section [T] this cannot be simply related to the 
mass and /or luminosity of the galaxies. So the trend shown 
in Figure [5] cannot simply be related a different spatial dis- 
tribution for g alaxies in differe nt luminosity bins, as done 
for example in lLin et alj (|2004l . see their figure 8). 



4 EVOLUTION OF SUBSTRUCTURES 

In this section, we study the evolution of substructures as a 
function of time, focusing in particular on their mass accre- 
tion histories and merger histories. In order to obtain these 
information, we have constructed merger histories for all 
self-bound haloes in our s i mulat ions, following the method 
adopted in ISpringel et al.1 (120051) and the improvements de- 
scribed in iDe Lucia fc BlaizotI l|2007[ ). 

Briefly, the merger tree is constructed by identifying 
a unique descendant for each substructure. For each sub- 
halo, we find all haloes that contain its particles in the 
following snapshot, and then count the particles by giving 
higher weight to those that are more tightly bound to the 
halo under consideration. The halo that contains the largest 
(weighted) number of its particles is selected as descendant. 
Next, all the pointers to the progenitors are constructed. By 
default, the most massive progenitor at each node of the 
tree is selected as the mam progenitor. iDe Lucia fc BlaizotI 
(|2007l ) noted that this can lead to ambiguous selections 
when, for example, there are two subhaloes of similar mass. 
In order to avoid occasional failures in the merger tree con- 
struction algorithm, they modified the definition of the main 
progenitor by selecting the branch that accounts for most of 
the mass of the final system, for the longest time. We have 
applied this modification to our merger trees. In this sec- 
tion, we consider only substructures that contain at least 
100 bound particles, and in a few cases, we use particular 
mass ranges to ease the comparison with the literature. 

In this section we will also study if the accretion and 
merger history of substructures depend on the environment, 
that we will approximate using the parent halo mass. It is 
worth stressing, however, that our haloes provide likely a 
biased sample for this analysis. In fact, excluding the most 
massive sample and some haloes that belong to the sample 
S2, all the other haloes reside in the regions surrounding the 
most massive haloes, which might not represent the 'typical' 
environment for halo in the same mass range. 



4.1 Mass Accretion History 

In this section, we use the merger trees constructed for 
our cluster sample to study the mass accretion histories 
of subhaloes of different mass and resi ding in di f ferent 
envir onments. Several previous studies (|Gao et al ] 12004 
iDe Lucia et al. 2004; Warnick et al. 2008) have pointed out 
that once haloes are accreted onto larger systems (i.e. they 
become substructures), their mass is significantly reduced 
by tidal stripping. The longer the substructure spends in a 
more massive halo, the larger is the destructive effect of tidal 
stripping. Previous studies have found that the efficiency of 
tidal stripping is large l y independent o f the parent halo mass 
l|De Lucia et al.ll2004l : Icao et al.|[20oi ) . 

We re-address these issues using all substructures re- 
siding within the virial radius of our haloes, and with mass 



larger than 10^'^ /i~^Mq at redshift z = (in our simula- 
tions, these substructures contain at least 100 particles). By 
walking their merger trees, following the main progenitor 
branch, we construct the mass accretion history (MAH) for 
all of these subhaloes, and record the accretion time {zaccr) 
as the last time the halo is a central halo, i.e. before it is ac- 
creted onto a larger structure and becomes a proper subhalo. 
Our final sample includes 39005 haloes, that we split in two 
bins of different mass by using either their present day mass 
or their mass at the accretion time. We end up with 33576 
haloes with mass larger than 10^^ h~^MQ at present (25344 
when using the mass at the accretion time) , and 5429 haloes 
with mass lower than the adopted threshold (13661 if the ac- 
cretion mass is used). In order to analyse the environmental 
dependence of the mass accretion history, we consider sep- 
arately subhaloes residing in our S5 and SI samples (these 
correspond to our lowest and largest parent halo mass, re- 
spectively). 



The top panels in Figure [U] show the distribution of the 
accretion times for the two mass bins considered. Left and 
right panels correspond to a splitting in mass done on the 
basis of the present day mass and of the mass at accre- 
tion, respectively. When considering the present day mass 
(left panel), the differences between the two distributions 
are small, with only a slightly lower fraction of more mas- 
sive substructures being accreted very late, and a slightly 
larger fraction of substructures in the same mass range be- 
ing accreted between z ~ 0.1 and z 1. A larger differ- 
ence between the two distribution can be seen when con- 
sidering the mass at the time of accretion (right panel). 
Substructures that are less massive at the time of accre- 
tion have been accreted on average later than their more 
massive counterparts. In particular, about 90 per cent of 
the substructures in the least massive bin considered have 
been accreted below redshift 0.5, while only 50 per cent of 
the most massive substructures have been accreted over the 
same redshift range. The distribution obtained for the most 
massive substructures is broader, extending up to redshift 
~ 2. This is largely a selection effect, due to the fact that 
we are only considering substructures that are still present 
at z = 0. Once accreted onto larger systems, substructures 
are strongly affected by tidal stripping so that, among those 
that were accreted at early times, only the most massive 
ones will still retain enough bound particles at present to 
enter our samples. The less massive substructures that were 
accreted at early times, have been stripped below the res- 
olution of our simulations and therefore do not show up in 
the solid histogram that is shown in the top right panel of 
Figure [S] 



The bottom panels of Figure |S] show the distribution of 
the ratios between present day mass and mass at accretion 
for subhaloes of different present day mass (left panel) and 
for different mass at accretion (right panel). Less massive 
subhaloes, which were accreted on average more recently, 
lose on average smaller fractions of their mass compared to 
more massive subhaloes for which the distribution is skewed 
to higher values. The difference between these distributions 
becomes more evident when one split the samples according 
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Figure 6. Left panels: distribution of the accretion times (top panel) and of the fraction of mass loss since accretion (bottom panel) 
for subhaloes of different mass at present time (different linestyles, as indicated in the legend). Right panels show the same distributions 
but for subhaloes split according to their mass at the accretion time. 
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Figure 7. Distribution of mass loss (ratio between the present 
day mass and the mass at accretion) for two different accretion 
ranges: solid line for Zaccr ^ 1 and dotted line for Zaccr < 1- 



to the mass at the time of accretion, as shown in the right 
paneL As explained above, however, this is affected by the 
fact that many of the least massive substructure will be 
stripped below the resolution of the simulation at z = 0. 
We have repeated the analysis done in Fig. [6] for subhaloes 



in each of the five samples used in our study, and we found 
there is no significant dependency on the environment. 

Fig. [7] shows that, as expected, substructures accreted 
earlier suffered significantly more stripping than substruc- 
tures that were accreted at later times. In particular, about 
90 per cent of subhaloes accreted at redshift larger than 1 
have been stripped by more than 80 per cent of their mass 
at accretion. For haloes that have been accreted at redshift 
lower than 1, the distribution is much broader, it peaks at 
~ 0.6 (i.e. about 40 per cent of the mass has been stripped 
for about 20 per cent of these haloes) but has a long tail to 
much lower values. Similarly to Fig. (6] we also tried to split 
this plot for different parent halo masses, without finding 
any significant trend with the environment. 



Fig. [8] shows the MAHs of subhaloes accreted at dif- 
ferent times. It shows results when subhaloes are split ac- 
cording to their present day mass (left panels), and the 
mass at accretion (right panels). As shown in previous stud- 
ies, the longer the halo is a substructure, the larger is its 
stripped mass. When substructures are split according to 
their present day mass, the influence of tidal stripping does 
not appear to depend strongly on the substructure mass. 
In contrast, if the mass at the accretion time is considered, 
in a given range of accretion times, haloes that are more 
massive lose a larger fraction of their mass with respect to 
their less massive counterparts. This is due to the fact that 
more massive haloes sink more rapidly towards the centre 
because of dynamical friction, and therefore suffer a more 
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significant stripping due to tidal interactions with the par- 
ent halo. Once again, this entails the fact that luminosity 
must correlate stronger with the subhalo mass computed at 
the time of accretion, i.e. before stripping had time to oper- 
ate. 



appear to be less stripped in the sample S5 than in SI (com- 
pare dashed and solid lines). The difference, however, is not 
large, but this might be affected by the fact that our haloes 
all reside in the regions surrounding very massive clusters. 



In Fig. |5] we plot the mean MAHs for subhaloes in 
the two mass bins considered and for two different 'envi- 
ronments', parametrized as the mass of the parent halo. In 
particular, we consider the samples S5 and SI (i.e. the least 
and the most massive haloes used in our study). Dashed 
and long-dashed lines show the MAHs for subhaloes in the 
sample S5 with mass in the range [lO'^" — 10^^] h~^MQ and 
larger than 10^^ h~^MQ, respectively. Solid and dotted lines 
show the MAH for subhaloes in the same mass ranges but 
for the sample SI. Here we consider the present day subhalo 
mass. Computing the same plot by adopting the subhalo 
mass at the time of accretion does not alter the results. We 
find that the environment does not significantly influence 
the mass accretion history of substructures. In the bottom 
panel, the long dashed line (corresponding to substructures 
more massive than 10^^ Mq in the sample S5) is likely 
affected by low number statistics. In the same panel a small 
difference can be seen for the less massive substructures that 



4.2 Merging Rate 

In recent years, a large body of observational evidence 
has been collected that demonstrates that galaxy interac- 
tions and mergers play an important role in galaxy evo- 
lution. In particular, numerical simulations have shown 
that major mergers between two spiral galaxies of com- 
parable mass can completely destroy the stellar disk and 
leave a kinematically hot remnant with structural and kine- 
matical properties similar to tho se of elliptical galaxies 
van den Bosch fc Whitell2010l . and references therein). 
Minor mergers and rapid repeated encounters wit h other 
galaxies residing in the same halo {harassment; Moore et al.l 
1 19961 : Im oore et akll 19981 ) can induce disk instabilities and /or 
the formation of a stellar bar, each of which affects the mor- 
phology of galaxies falling onto clusters. As galaxy mergers 
are driven by mergers of the parent dark matter haloes, it 
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Figure 9. Average mass accretion history for subhaloes in three 
different ranges of accretion time, as a function of environment. 
Dashed and long-dashed hnes show the MAH for subhaloes in the 
sample S5 with mass in the range [10^" — 10^^] h~^MQ and larger 
than 10^^ h~^MQ, respectively. Solid and dotted lines show the 
MAH for subhaloes in the same mass ranges but for the sample 
SI. 



is interesting to analyse in more detail the merger statistics 
of dark matter substructures. 

The mass accretion history discussed in the previous 
section does not distinguish between merger events (of dif- 
ferent mass ratios) and accretion of 'diffuse material'. In 
order to address this issue, and in particular to study the 
merger rates of dark matter substructures, we have taken 
advantage of the merger trees constructed for our sam- 
ples. We have selected all subhaloes with mass larger than 
10^^ h-^MQ at redshift zero, and have followed them back 
in time by tracing their main progenitor branch, and record- 
ing all merger events with other structures. In particular, we 
take into account only mergers with objects of mass larger 
than 10^° h~^MQ, and mass ratios larger than 5 : 1. We note 
that both these values are computed at the time the halo is 
for the last time central (the mass of the main progenitor 
at the time of accretion is considered to compute the mass 
ratio). 

Fig. [To] shows the merging rate for all subhaloes that 
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Figure 10. Mean number of major mergers as a function of red- 
shift, for subhaloes in two different ranges of accretion time. We 
take into account only subhaloes with mass M ^ IQi^ h-'^MQ at 
redshift z = and merger events that include systems with mass 
M > IQio /i-iM0. 
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Figure 11. Mean number of major mergers as a function of red- 
shift, for subhaloes in different environments, quantified as the 
mass of their parent halo. As in Fig. llOl we take into account only 
subhaloes with mass M ^ IQI^/j-IMq at redshift z = Q and 
merger events that include systems with mass M ^ 10^'' /i~^Mq. 
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satisfy the above conditions. We consider in this plot only 
objects that experienced at least one merger event. The solid 
line shows the mean number of mergers for subhaloes that 
were accreted at z < 0.5, while the dotted line shows the 
resulting merger rate for objects accreted between 0.5 < z ^ 
1. The figure shows that in both cases, the slope of the lines 
become shallower close to the accretion time, i.e. mergers 
between substructures are suppressed because of the large 
velocity dispersion of the parent haloes. Interestingly, haloes 
that were accreted earlier experience, on average, one more 
major merger than haloes accreted at later times. 

We repeat the same analysis looking at the merging rate 
as a function of environment. Fig. [TT] shows the cumulative 
number of mergers for subhaloes in our five samples. The 
mean number of mergers increases as a function of the parent 
halo mass, although subhaloes in the sample S4 experience 
on average fewer mergers than subhaloes in the sample S5. 
This is not surprising since subhaloes in the surroundings of 
more massive haloes have a larger probability to merge with 
other structures. 



5 DISCUSSION AND CONCLUSIONS 

We have used a large set of high-resolution simulated haloes 
to analyse the statistics of subhaloes in dark matter haloes, 
and their dependency as a function of the parent halo mass 
and physical properties of the parent halo. While some of 
the results discussed in this study confirm results from pre- 
vious studies, it is the first time that a systematic analysis 
of the properties and evolution of dark matter substructures 
is carried out using a large simulation set carried out using 
the same cosmological parameters and simulation code. Our 
main results can be summarized as follows: 

(i) More massive haloes contain increasing fractions of 
mass in subhaloes. This does not exceed ~ 10 per cent of 
the total mass, in agreement with previous studies. There is, 
however, a very large halo-to-halo scatter that can be par- 
tially explained by a range of halo physical properties, e.g. 
the concentration. Indeed, in more concentrated haloes sub- 
structures suffer of a stronger tidal stripping so that they are 
characterized by lower fractions of mass in substructures. 

(ii) We find that the subhalo mass function depends weakly 
on the parent halo mass and on redshift. This can be ex- 
plained by considering that haloes of larger mass are less 
concentrated and dynamically younger than their less mas- 
sive counterparts, and that haloes of a given mass are on 
average less concentrated at higher redshift. Our fin dings 
confirm results from previous studies (|Gao et al.ll201ll ). and 
extend them to larger halo masses. 

(iii) As shown in p revious work (e.g. iGhigna et all 1 19981 : 
iDe Lucia et al.ir2004 ) , subhaloes are anti-biased with respect 
to the dark matter in the inner regions of haloes. The anti- 
bias is considerably reduced once subhaloes are selected on 
the basis of their mass at the time of accretion, or neglect- 
ing those that were accreted at later times. We also find 
that the spatial distribution of subhaloes does not depend 
significantly on halo m ass, as suggested in previous work by 
iDe Lucia et al.l (|2004l ). The most massive substructures are 
located at the outskirts of haloes and this mass segregation 
is more important at higher redshift. 



(iv) Once accreted onto larger systems, haloes are strongly 
affected by tidal stripping. The strength of this stripping ap- 
pears to depend on the mass of the accreting substructures: 
those that are more massive at the time of accretion tend to 
be stripped by larger fractions of their initial mass. 

(v) Mergers between substructures are rare events. Follow- 
ing the merger trees of substructures, however, we find that 
they have suffered in the past about 4-5 important (mass 
ratio 1:5) mergers. As expected, the number of mergers ex- 
perienced depends on the environment: subhaloes in more 
massive systems have experienced more mergers than those 
of similar mass residing in less massive haloes. 

Dark matter substructures mark the sites where lumi- 
nous satellites are expected to be found, so their evolution 
and properties do provide important information on the 
galaxy population that forms in hierarchical models. As dis- 
cussed in previous studies, however, because of the strong 
tidal stripping suffered by haloes falling onto larger struc- 
tures, it is not possible to simply correlate the population 
of subhaloes identified at a given cosmic epoch to that of 
the corresponding galaxies. The galaxy luminosity/stellar 
mass is expected to be more strongly related to the mass 
of the substructure at the time of infall and, depending on 
the resolution of the simulations, there might be a significant 
fraction of the galaxy population that cannot be traced with 
dark matter substructures because they have been stripped 
below the resolution limi t of the sim ulation (the 'orphan' 
galaxies - see for example IWang et al.| [2006). 

Nevertheless, our results do provide indications about 
the properties of the galaxy populations predicted by hierar- 
chical models. Tidal stripping is largely independent of the 
environment (we have parametrized this as the parent halo 
mass), while the accretion rates of new subhaloes increases 
at increasing redshift. The nearly invariance of the sub- 
halo mass function results from the balance between these 
two physical processes. If the amount of dark matter sub- 
structures is tracing the fraction of recently infallen galax- 
ies, the fraction of star forming galaxies is expected to in- 
crease with increasing redshift (the 'Butch er-Oemler' effect, 
iButcher fc Oemlerl [19781 , lKauffmann|[l995l ) . In addition, our 
findings suggest that stronger mass segregation should be 
found with increasing redshift. 

There is a large halo-to-halo scatter that can be only 
partially explained by a wide range of physical properties. 
This is expected to translate into a large scatter in e.g. the 
fraction of passive galaxies for haloes of the same mass, 
with more concentrated haloes hosting larger fraction of 
red/passive galaxies. Finally, there is an obvious merger bias 
that is expected to translate into a different morphological 
mix for haloes of different mass. In future work, we plan to 
carry out a more direct comparison with observational data 
at different cosmic times, by applying detailed semi-analytic 
model to the merger trees extracted from our simulations. 
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